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Abstract. A Brownian time process is a Markov process subordinated to the 
absolute value of an independent one-dimensional Brownian motion. Its transition 
densities solve an initial value problem involving the square of the generator of the 
original Markov process. An apparently unrelated class of processes, emerging as 
the scaling limits of continuous time random walks, involve subordination to the 
inverse or hitting time process of a classical stable subordinator. The resulting 
densities solve fractional Cauchy problems, an extension that involves fractional 
derivatives in time. In this paper, we will show a close and unexpected connection 
between these two classes of processes, and consequently, an equivalence between 
these two families of partial differential equations. 



1. Introduction 

The goal of this paper is to estabhsh a connection between two seemingly disparate 
classes of subordinated stochastic processes, and their governing equations. The first 
class of processes uses Brownian motion as a subordinator [U El UHl EHl EDI l2T] . 
The second class uses the inverse or hitting time process of a nondecreasing stable 
Levy process as a subordinator [38l [391 SOl [IS, HI] . One application of the Brownian 
subordinator involves scaling limits of transverse diffusion in a crack [2l]. Inverse 
Levy subordinators occur in scaling limits of continuous time random walks used to 
model anomalous diffusion in physics [l3l [Hj, [6l], finance [271 Il2l [57], and hydrol- 
ogy [ISl [151 159] • Neither class of processes is Markovian. However, their transition 
densities satisfy certain abstract differential equations that characterize the evolution 
of the process. Given a Markov process whose transition densities solve a Cauchy 
problem, a Brownian subordinator yields another process whose one dimensional dis- 
tributions solve a related initial value problem involving the square of the generator 
of the Markov semigroup. Subordinating the same original Markov process to an 
inverse stable Levy process leads to a fractional Cauchy problem, where the integer 
time derivative is replaced by a fractional derivative whose order equals the stable in- 
dex. Both situations lead to anomalous subdiffusion, where probability mass spreads 
slower than the classical t^^"^ rate seen in Brownian motion. In this paper we will 
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show that, when the order of the fractional derivative equals 1/2, the two govern- 
ing equations have the same unique solution, and hence the two processes have the 
same one dimensional distributions. Similar observations apply for other values of 
the fractional time derivative, and there are connections to subordination by a sym- 
metric stable process. In applications to physics, the principal appeal of both classes 
is that the subordinators have a pleasant scaling property, that can be inherited by 
the subordinated processes. 

A Levy process B{t) has stationary, independent increments and t ^ B{t) is con- 
tinuous in probability, see for example Bertoin [TB] or Sato [5Bj. Classical Brownian 
motion in is the special case where B{t) has a mean zero normal distribution. In 
this case, the probability distribution p{x,t) of x = B{t) solves the diffusion equation 
dp/dt = D Ap for some D > 0, where A = ^^d"^ / dxj. The random walk model of 
particle motion, where particles undergo independent identically distributed jumps at 
regularly spaced interval in time, converges to Brownian motion as both the spatial 
and time scales increase to infinity. This close connection between Brownian motion 
and the diffusion equation is due to Bachelier |5] and Einstein [23]. Sokolov and 
Klafter [53] discuss modern extensions to include heavy tailed particle jumps and 
random waiting times, leading to fractional diffusion equations, the simplest of which 
is d^p/dt^ = — D (— A)"/^p. The fractional Laplacian reflects heavy tailed particle 
jumps, where the probability of a jump magnitude exceeding some large r > di- 
minishes like r~", and the fractional time derivative codes heavy tailed waiting times, 
where the probability of waiting longer than t falls off like t~^. Heavy tails in space 
lead to super diffusion, where a plume of particles spreads faster than the classical 
t^/^ rate associated with Brownian motion. Heavy tails in time lead to subdiffusion, 
since long waiting times retard particle motion. See Metzler and Klafter [ISj HI] for 
a recent survey. 

A completely different model of particle motion emerges from the theory of random 
walks in random media, pioneered by Sinai [62] . These models lead to interesting lo- 
calization phenomena [25l [53] not seen in the case of random waiting times [17] . A 
closely related problem is diffusion on fractals [3 El Ej- One particular application, 
diffusion in a crack, led Burdzy and Khoshnevisan [2T] to consider Brownian subordi- 
nators. In particular, a two-sided Brownian motion subordinated to an independent 
Brownian subordinator yields the scaling limit of transverse diffusion in a crack. Es- 
sentially, the crack is modeled using the sample path of the outer process, and the 
inner process or subordinator codes the subdiffusive effect of restricted motion in the 
crack. Given this history, it would be difficult to imagine that these two classes of 
non-Markovian processes exhibit any close connections. However, we will show that 
they are very closely related. 
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2. Background 



Let B{t) be a Brownian motion on M!^ and Yt another independent one-dimensional 
Brownian motion. Allouba and Zheng [H |2] introduce a process they call Brownian 
time Brownian motion (BTBM) defined a.s Zt = B{\Yt\). A related process called iter- 
ated Brownian motion process (IBM) was first considered by Burdzy [18]. Take -Bi(t), 
B2{t) to be independent Brownian motions on R'^, and Yt as before. Define a two-sided 
Brownian motion by B{t) = Bi{t) for t > and B2{-t) for t < 0. Then the IBM 
process is defined by Zt = B{Yt). Various extensions of the BTBM or IBM have been 
considered, including a general Markov outer process [H |2] and a symmetric stable 
subordinator studied by Nane [ITJIIS]- Similarly, the excursion-based Brownian-time 
process is defined by breaking up the path of into excursion intervals (i.e., 

maximal intervals of time on which |V(t)| > 0) and, on each such interval, choosing 
an independent copy of the Markov process from a finite or an infinite collection. A 
simple conditioning argument shows that all these Brownian time processes have the 
same transition density functions, and hence the same one-dimensional distributions 
[2j . A similar remark holds for symmetric stable subordinators. That these processes 
are non-Markovian can be seen by noting that the transition densities do not solve 
the Kolmogorov equations. 

Allouba and Zheng [2J and DeBlassie [23] show that for iterated Brownian motion 
(IBM) Zt = BiYt) the function 

u{t,x) = E,[f{Zt)] := E[f{Zt)\Zo = x] 

solves the initial value problem 

(2.1) ^u{t,x) = ^^ + AMt,xy, u{0,x) = f{x) 



for t > and x G M"'. The non-Markovian property of IBM is refiected in the 
appearance of the time-variable initial term in the PDE. 

For a Markov process X, the family of linear operators T{t)f{x) = Ex[f{X{t))] = 
E[f{X{t))\X{0) = x] forms a bounded continuous semigroup on the Banach space 
L^(M'^), and the generator L^f^x) = limhiQ h~^{T{h)f{x) — f{x)) is defined on a 
dense subset of that space, see for example [H [2H]- Then u{t,x) = T(t)f{x) solves 
the Cauchy problem 

(2.2) ^u{t, x) = LJix); u{0, x) = f{x) 

for t > and x G M'^. Allouba and Zheng [2] show that if we replace the outer process 
B{t) with a continuous Markov process X{t), the same result holds, except that 
we replace the Laplacian in the PDE (12.11) with the generator of the continuous 
semigroup associated with this Markov process. That is, u{t,x) = Ex[f{Zt)] solves 
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the initial value problem 

(2.3) ^u{t, x) = h^il^ + L,\{t, x)- m(0, x) = fix) 
for t > and x G M'^. 

Remark 2.1. Iterated Brownian motion (IBM) has many properties analogous to 
those of Brownian motion. The process Zt has stationary (but not independent) 
increments, and is self-similar with index 1/4, meaning that Zct and (^I'^Zt have the 
same finite dimensional distributions for every c > 0. Burdzy (1993) [18] showed that 
IBM satisfies a Law of the iterated logarithm (LIL): 

Zt 2^/4 
limsup —rTTT, — ; — -, — , ^^o/. = -^m cl-s. 

t^oo^tl/^(loglog(l/t))3/4 33/4 

A Chung-type LIL by was established by Khoshnevisan and Lewis [35j and Hu et 
al. [30]. Khoshnevisan and Lewis |36] extended results of Burdzy [T9| to develop a 
stochastic calculus for IBM. Local times of IBM were studied by Burdzy and Khos- 
nevisan [20j, Csaki et al. [22], Shi and Yor [60j, Xiao [M], and Hu [2H]. Banuelos and 
DeBlassie [TT] studied the distribution of exit place for IBM in cones. Nane studied 
the lifetime asymptotics of IBM on bounded and unbounded domains [121 SHI 150] . 
and generalized isoperimetric-type inequalities to IBM [49j . 

Remark 2.2. Funaki [26j considered a different version of iterated Brownian motion, 
and studied the PDE connection of that process. Nane [17] established the PDE 
connection of a related class of processes defined in a similar manner, replacing the 
inner time process with a symmetric stable process. 

Zaslavsky [65] introduced the fractional kinetic equation 

(2.4) ■^u{t,x) = L^u{t,x); u{0,x) = f{x) 

for Hamiltonian chaos, where < /9 < 1 and is the generator of some continuous 
Markov process Xo(t) started at x = 0. Here d^g(t)/dt^ is the Caputo fractional 
derivative in time, which can be defined as the inverse Laplace transform of s^g{s) — 
s^~^g{0), with g{s) = e~^*g{t)dt the usual Laplace transform. Baeumer and 
Meerschaert [TD] and Meerschaert and SchefHer [ID] show that the fractional Cauchy 
problem (12. 4p is related to a certain class of subordinated stochastic processes. Take 
Dt to be the stable subordinator, a Levy process with strictly increasing sample paths 
such that E[e~^^^] = e"**'', see for example [HI [56]. Define the inverse or hitting time 
or first passage time process 

(2.5) Et = inf{x > : D{x) > t}. 

The subordinated process Zt = Xo^Et) occurs as the scaling limit of a continuous 
time random walk (also called a renewal reward process), in which iid random jumps 



are separated by iid positive waiting times [ID]. If the waiting times Jj satisfy P{Ji > 
t) = t~^L(t) where < /3 < 1 and L{t) is slowly varying, then they belong to the 
strict domain of attraction of a stable law with index /?, and their partial sum process 
St = Ji + ■ ■ ■ + J[t] converges after rescaling to the process Dt, in the Skorokhod Ji 
topology on D(]R+,R). The number of jumps by time t > is given by the inverse 
process Nt = max{n : S'„ < t} and |1D] shows that a rescaled version of Nt converges 
to the hitting time process Et in the same topology. Similarly, if the iid vector particle 
jumps Yi independent of the waiting times satisfy a multivariable regular variation 
condition (or if they have a finite covariance matrix), then the partial sum process 
V{t) = Yi + ■ ■ ■ + Y[t] converges after linear operator rescaling to the operator Levy 
motion Xo(t), in the Skorokhod Ji topology on D(R+,R'^), see |10l Theorem 4.1]. 
An operator Levy motion is a Levy process such that Xo(t) has a centered operator 
stable distribution, see Example 11.2.18]. This means that Xo(ct) and c^Xo(t) 
are identically distributed for all c > 0, for some linear operator E, see for example 
[SUEZ]. The continuous time random walk V{St) models the location of a particle at 
time t > 0. A continuous mapping argument, under some mild technical conditions, 
yields convergence of the rescaled CTRW to the subordinated process Xo{Et) in 
the somewhat weaker Mi topology, see ^40]. Given a Banach space and a bounded 
continuous semigroup T{t) on that space with generator L^, it is well known that 
p(t,x) = T(t)f{x) is the unique solution to the abstract Cauchy problem 



for any / in the domain of see for example [U [52]. Theorem 3.1 in shows that, 
in this case, the formula 



yields the unique strong solution to the fractional Cauchy problem (12.41) . Here (7/3 (t) 
is the smooth density of the stable subordinator, such that the Laplace transform 
gp[s) = /~ e-'^gp{t) dt = e-'^. Choose x eW^ and let X{t) =x + Xo{t). It follows 
from Theorem 5.1 in pU] that, in the special case where T{t)f{x) = Ex[f{X{t))] is 
the semigroup on Li(R'^) associated with the operator Levy motion X{t), the same 
formula (12. 7p also equals u(t,x) = Ex[f{Zt)] where Zt = X{Et) is the CTRW scaling 
limit process. Hence the subordinated process Zf is the stochastic solution to the 
fractional Cauchy problem (12. 4p in this case. 

In the case where Xo(t) is a Levy process started at zero and X(t) = x + Xo(t) 
for X G M'^, the generator L^. of the semigroup T(t)f{x) = Ej.[f{X(t))] is a pseudo- 
differential operator [3l [3T[ [58] that can be explicitly computed by inverting the Levy 
representation. The Levy process Xo(t) has characteristic function 



(2.6) 



d_ 
di 



■pit, x) = L^p{t, x); p(0, x) = f{x) 



(2.7) 




E[exp{ik ■ Xo{t))] = exp{tij{k)) 
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with 

^(A;) =^k-a-h-Qk+ [ (e^^-y - 1 - ^ 

where a G M'', Q is a nonnegative definite matrix, and is a cr-finite Borel measure 
on such that 

/ min{l, < oo, 

see for example 071 Theorem 3.1.11] and [S] Theorem 1.2.14]. Let 

denote the Fourier transform. Theorem 3.1 in [10] shows that Lxf{x) is the inverse 
Fourier transform of ip{k)f{k) for all / G D{Lj.), where 

D{L,) = {/ G L\R'') : ^(A;)/>) = h{k) 3 he L\R'')}, 
LJ{x) = a ■ V/(x) + Iv ■ QV/(x) 

+ / f /(^ + - / w - n^^) m) 

Jy^o \ i + y^ J 

for all / G iy^'^(]R'^), the Sobolev space of L^-functions whose first and second partial 
derivatives are all L^-functions. This includes the special case where Xo{t) is an oper- 
ator Levy motion. We can also write = ip{—iS/) where V = {d/dxi, . . . , d/dxd)' ■ 
For example, if Xo(t) is spherically symmetric stable then ip{k) = —D\\k\\°' and 
Lx = — D(— A)°/^, a fractional derivative in space, using the correspondence kj 
—id/dxj for I < j < d. If Xq has independent stable marginals, then one possible 
form is ip{k) = DY2j{ikj)°'^ and = D'Y^.d°'^/dx'^ using Riemann-Liouville frac- 
tional derivatives in each variable. This form does not coincide with the fractional 
Laplacian unless all aj = 2. 

Remark 2.3. The literature on fractional calculus uses a different semigroup definition. 
For example, in piTj the semigroup associated with a Levy process Xo(t) started at 
x = is defined by 

TFc{t)f{x) = E[f{x - Xoim = I fi^- y)Pxo(t){dy). 

One physical interpretation of this formula, when f{x) is the probability density of a 
random variable W, is that W represents the location of a randomly selected particle 
at time t = 0, and Tpc{t)f{x) is the probability density of the random particle 
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and 
(2.8) 



location W + XQ{t) at time t > 0. In this paper, we use the semigroup definition from 
the Markov process hterature, based on the process X(t) = x + Xo(t), so that 



T{t)f{x) = E4f{X{t))] = / f{x + y)Px,it){dy). 

Clearly this is just a matter of replacing Xq by — Xq. The paper [TO] also uses a 
different definition fpcik) = J e^^'^ f{x) dx for the Fourier transform. Each of these 
two changes implies a change of A; to —/c in the formula for the Fourier transform of 
the semigroup, and hence the Fourier symbol '^(fc) is the same for both. However, 
the interpretation of the Fourier symbol as a pseudo-differential operator changes. In 
the notation of this paper, the derivative V/(x) has Fourier transform {ik)f{k) but 
in the notation of [10] the Fourier transform is {—ik)f{k). Hence the generator = 
ipliV) in that paper, which explains why equation (12.81) differs from the corresponding 
formula (8) in [TO] . 

Any Markov process semigroup operator T(t)f{x) = Ex[f{X(t))] is a pseudo- 
differential operator 

T(t)/(x) = (27r)-<^ / e''-^\t{x,k)f{k)dk 



on the space of rapidly decreasing functions [33. Theorem 1.4]. Under some mild 
conditions, if the domain D{Lx) of the generator of the extension of this semigroup 
to the Banach space of bounded continuous functions contains all smooth functions 
with compact support, then one can write 

LJ{x) = a{x) ■ V/(x) + • Q{x)Vf{x) 

^^'^^ + / (fix + y)- fix) - ^P^) x) 

where a(x) G M'^, Qix) is a nonnegative definite matrix, and (j)idy,x) is a a-finite 
Borel measure on for each x G M'^ such that 



/ mm{l, \ \y\\'^}(j){dy, x) < oo, 



see [32]. In this case, the Fourier symbol 

^(fc, x) = zk- aix) - h ■ Q{x)k + [ (e'^ y - 1 - 4>{dy, x), 

and Lj.f{x) = ■?/'(— iV, is the inverse Fourier transform of 4'ik,x)fik). 

Remark 2.4. A connection between time-fractional equations and Brownian time was 
also noticed by Orsingher and Behgin [51j , using a very different approach. They show 
that the density functions of iterated Brownian motion solve the one-dimensional 
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time-fractional equation d^^'^ujdi}'!'^ = d'^u/dx^ by considering this equation as the 
end-member of a class of fractional telegraph equations. 



3. Main results 

The results of this section establish a connection between two seemingly distinct 
classes of subordinated stochastic processes. Markov processes are stochastic solu- 
tions to the abstract Cauchy problem fl2.6p . Brownian subordinators yield stochastic 
solutions to an initial value problem (12.31) involving the square of the Markov genera- 
tor. Inverse stable subordinators yield solutions to a fractional Cauchy problem (12.41) 
with the same Markov generator. In the case where the stable index (3 = 1/2, our 
first results shows that these two equations have the same solutions, and hence, the 
subordinated processes have the same one- dimensional distributions. 

Theorem 3.1. Let be the generator of a continuous Markov semigroupT{t)f{x) = 
Ex[f{Xt)], and take f G D{Lx) the domain of the generator. Then, both the initial 
value problem (12.31) . and the fractional Cauchy problem (12. 4p with (3=1/2, have the 
same solution given by 

(3.1) u{t,x) = -|=^"T(.)/(x)exp (^-^^ ds. 

Proof. Note that the Markov semigroup T{t)f{x) = Ex[f{Xt)] is a uniformly bounded 
and strongly continuous semigroup on the Banach space L^(M'^). Then Theorem 3.1 
in Baeumer and Meerschaert \10\ shows that the function u{t,x) defined by (12.71) is 
the unique solution of the fractional Cauchy problem (12. 4p for any / G D[Lx). In the 
case P = 1/2, we have Example 1.3.19] 

1 /I 



and then a change of variables shows that 

t t f°° 

(3.2) u{t,x) = - p{x,s)gp{-^)s'^"^~'^ds= I p{x, s)q{t, s)ds 

where p{x,t) = T{t)f{x) and 

2t ^ s'^\ 



(3.3) 



'v/4^ 
: exp 



'At 



is a probability density on s > for all t > 0. Hence we have that (13. ip is the unique 
solution to the fractional Cauchy problem (12.40 . 

Allouba and Zheng [2j show that the initial value problem (12. 3 p has solution 
u{t, x) = Ex[f{Zt)\ where Zt = B{\Yt\) is the IBM process. It is not hard to check that 



the function q{t, s) in fl3.3p is the probabihty density of the Brownian subordinator 
\Yt\. Then a simple conditioning argument shows that 

POO 

(3.4) uit,x) = E4f{Zt)]= qit,s)pix,s)ds 

Jo 

where p{x,s) = T(s)/(x) = E^f^Xf)] is the unique solution to the initial value 
problem (12.61) . Hence both the fractional Cauchy problem (12. 4p and the initial value 
problem (12.31) have the same solution. □ 

Corollary 3.2. For f G D(A^), both the initial value problem 

(3.5) ^u{t,x) = ^^ + Aluit,xy, uiO,x) = fix) 



and the fractional Cauchy problem 

51/2 

(3.6) -^^uit^x) = A^u{t,x)] u{0,x) = f{x) 
have the same solution given by (12.71) . where 

(3.7) p{t,x) = T{t)f{x) = J^^ fix - y){A7Tt)-''/' exp dy. 

Proof. Apply Theorem 13.11 and note that (13.71) is the density of the Brownian motion 
Markov process with generator A^. □ 

Suppose that X(t) is an operator Levy motion, Ft = inf{x > : > t} is the 
inverse or hitting time process of the stable subordinator Dt with F[e~^^*] = e"*'^'', 
and Zf = X{Ft). Let F^ be the generator of the semigroup T(t)f{x) = F4f{X(t))], 
a pseudo-differential operator defined by (12.81) . It follows from Theorem 5.1 in |40] 
that, for any initial condition / G D{Fx), the function u(t,x) = F^f^Zt)] solves 
the fractional Cauchy problem (12.41) . Hence we call the non-Markovian process Zt 
the stochastic solution to this abstract partial differential equation. The following 
theorem extends this result to a broader class of driving processes. 

Theorem 3.3. For any continuous Markov semigroup X{t) with generator F^, let 
Ft = ini{x > : > t} be the inverse or hitting time process of the stable 
subordinator Dt, independent of X , with £'[e~'^^*] = e"*^'' for some < P < 1, 
and take Zt = X{Ft). Then for any initial condition f G Di^F^), the function 
u{t,x) = Fj.[f{Zt)] solves the fractional Cauchy problem (12.41) . 

Proof. The proof is very similar to Theorem 13.11 above. Theorem 3.1 in [10] shows 
that the function u{t,x) defined by (12.71) is the (unique) solution of the fractional 
Cauchy problem (12.41) for any / G D{Fx). A change of variables shows that (13.21) 
holds, where gp is the density of the random variable Di. Corollary 3.1 in [40j shows 
that the function q{t, s) = tl3~^gi3{s~^/^t)s~^/^~^ is the density of the hitting time Ft., 
and then the result follows by a simple conditioning argument. □ 
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Corollary 3.4. For any continuous Markov process X[t), both the Brownian-time 
subordinated process X{\Yt\) and the process X{Et) subordinated to the inverse 1/2- 
stable subordinator have the same one- dimensional distributions. Hence they are both 
stochastic solutions to the fractional Cauchy problem (12.41) . or equivalently, to the 
higher order initial value problem (12.31) . 

Proof. In the proof of Theorem 13.11 we noted that the function s) in (13.31) is the 
probability density of the Brownian subordinator s = \Yt\ where Yt is a standard 
Brownian motion. The proof of Theorem 13.31 shows that this function is also the 
probability density of the inverse 1/2-stable subordinator s = Ef. Then the result 
follows from a simple conditioning argument. □ 

Remark 3.5. The proof of Theorem 13.11 relies on Theorem 0.1 in Allouba and Zheng 
[2], but the proof of Theorem 0.1 in that paper may not be completely rigorous. The 
essential argument, using integration by parts twice, is that 



which shows that (13.41) solves the higher order initial value problem (12.31) . In a later 
paper pj, Allouba points out that the last step where the operator is passed 
outside the integral is not obvious, and adds this as a technical condition, which is 
then verified in the special case = the Laplacian operator. 

The next result is a restatement of Theorem 13.11 for Levy semigroups. The proof 
does not use Theorem 0.1 in Allouba and Zheng p], rather it relies on a Laplace- 
Fourier transform argument. We will use the following notation for the Laplace, 
Fourier, and Fourier-Laplace transforms (respectively): 



d_ 
dt 



u(t, x) 






e ^^u{t, x)dtdx. 
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Theorem 3.6. Suppose that X{t) = a; + Xo(t) where XQ{t) is a Levy process starting 
at zero. If L^. is the generator fl2.8l) of the semigroup T{t)f{x) = on 
L^{W^), then for any f G D{L^), both the initial value problem (12.31) . and the frac- 
tional Cauchy problem (12.41) with (3 = 1/2, have the same unique solution given by 
(O). 



Proof. Take Fourier transforms on both sides of (12.31) to get 
(3.8) = -l=i;(k)f{k) + m'Ht, k) 



using the fact that ip{k)f{k) is the Fourier transform of L^f^x). Then take Laplace 
transforms on both sides to get 

su{s, k) - u{t = 0, A;) = s-^/^tlj{k)f{k) + i){kfu{s, k), 

using the well-known Laplace transform formula 

f-l3 



r(i-/3f * = ^ 

for P < 1. Since u(t = 0,k) = f{k), collecting like terms yields 

,3.9) (i±£^^g)M 

for s > sufficiently large. 

On the other hand, taking Fourier transforms on both sides of (12.41) with (3 = 1/2 
gives 

(3.10) ^^^^ = mu{t,k) 

Take Laplace transforms on both sides, using the fact that s^g{s) — s^^^g{0) is the 
Laplace transform of the Caputo fractional derivative d^^ g{t) / dt^ , to get 

s^'^u{s, k) - s-^'^f{k) = ^{k)u{s, k) 
and collect terms to obtain 

s-'l^f\k) 



u{s, k) 
(3.11) 



SV2 _ ^(A;) 

s-^'^f\k) s^l^ + ^{k) 

(l + ,-l/2^(fc))/>) 



s — ip{ky 

which agrees with (13.91) . For any fixed k G W^, the two formulas are well-defined and 
equal for all s > sufficiently large. 
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Theorem 3.1 in Baeumer and Meerschaert [TU] implies that the function u{t,x) 
defined by (12.71) solves the fractional Cauchy problem (12.41) in L^(M'^) for any / G 
D{Lx). Take Fourier transforms in (12. 7p and apply the Fubini theorem to get 

POO 

(3.12) n(t, k) = / />) exp((t/s)^/2^(A;))(?i/2(s) ds. 

Jo 

Note that exp{{t/ sy^'^%jj{k)) is bounded since ijjik) is negative definite, and then a 
simple dominated convergence argument shows that u{t, k) is continuous in t > for 
any k G M*^. Hence the uniqueness theorem for Laplace transforms ^ Theorem 1.7.3] 
shows that, for each k G W^, (I3.12p is the unique continuous function whose Laplace 
transform is given by either (13. 9p or (13. lip . Since x ^ u{t, x) is an element of L^(R°') 
for every t > 0, and since two elements of L^(R'^) with the same Fourier transform 
are equal cix-almost everywhere, (13. ip is the unique element of L^(M'^) whose Fourier 
transform equals (I3.12p . 

Now if u{t, x) is any solution to (12. 3p . then it has Fourier-Laplace transform u{s, k) 
given by (13.90 or equivalently by (I3.1ip . Since u(t,k) solves (13. Sp . it is differentiable 
in t > and hence continuous. Then the above argument shows that this solution 
is equal to (12. 7p rfx-almost everywhere for every t > 0. Hence (12.70 solves both the 
higher order Cauchy problem (12.30 . and the fractional Cauchy problem (12. 4p with 
(3 = 1/2, for any / G D{L^), and it is the unique solution as an element of L^(]R'^), 
i.e., it is unique a.e.-dx for every t > 0. The proof of Theorem 13.11 shows that (12.70 
reduces to (13. ip in this case. □ 

Remark 3.7. Uniqueness to (12. 4p was shown by Bajlekova in [6, Corollary 3.2]; Allouba 
and Zheng [2] did not discuss uniqueness of solutions. 

Remark 3.8. It is reasonable to conjecture that Theorem 13.61 can be extended to a 
sub-class of Markov processes X(t) whose generators are pseudo-differential operators 
with negative definite symbols ilj{k,x). One technical difficulty is that il){k,x)^ is 
usually not the Fourier symbol of a semigroup generator for integers n> 1. 

The Fourier-Laplace transform method can be extended to establish connections 
between other fractional Cauchy problems and their corresponding higher-order initial 
value problems. The next result gives one such correspondence. 

Theorem 3.9. Suppose that X{t) = x + XQ{t) where Xo(t) is a Levy process starting 
at zero. If is the generator (12. 8p of the semigroup T{t)f{x) = Ex[{f{Xt))], then 
for any f G D{Lx), both the initial value problem 

dv(t r) /~2/3 ^-1/3 

(3.13) ^^^^ = f^)LJ{x) + ^^^Llf{x) + Llu{t,xy, u{0,x) = f{x) 

and the fractional Cauchy problem (12.40 with f3 = 1/3, have the same unique solution 
given by (12. 7p where gi/3{t) is the probability density of the 1/3-stable subordinator, 

so that Jg°° e~'^^gi/s{t) dt = e~'^^^^ for all s > 0. 
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Proof. The proof is very similar to Theorem 13.61 Take Fourier transforms on both 
sides of f l3.13p to get 

= ^^^(^)/(^) + f(273)^^^^'^^^^ + 
then take Laplace transforms on both sides to get 

su{s, k) - f{k) = s-^'^ij{k)f{k) + s-^'^4j{kff{k) + ij{kfu{s, k) 
and collect terms to obtain 

(3.14) ^^^^^^+s-y^m^s^^^i>ef)m 

for s > sufficiently large. 

On the other hand, taking Fourier transforms on both sides of (12.41) with (3=1/3 
gives 

■^u{t,k) = ij{k)u{t,k), 
and then taking Laplace transforms yields 

s^/^u{s, k) - s~^l^j(k) = ^{k)u{s, k). 

Collecting terms yields 

s~ 



u{s, k) 
(3.15) 



s^l^-i){k) 



s^^-tpik) s2/3 + si/3^(A;) + ^(A;)2 
(1 + s-i/3^(A;) + s-2/3^(A;)2)/(A;) 



s — Tpiky 

which agrees with (13.141) . For any fixed k G M'^, the two formulas are well-defined 
and equal for all s > sufficiently large. The remainder of the argument is identical 
to Theorem 13.61 □ 

Corollary 3.10. Suppose thatX{t) = x+Xo{t) where XQ{t) is a Levy process starting 
at zero, and let Et = inf {a; > : Dx > t} be the inverse or hitting time process of the 
1/3-stable subordinator Dt, independent of X , with E[e~^^^] = e~^^^'^ . If is the 
generator (12.81) of the semigroup T(t)f{x) = E^Kf^Xt))] and Zt = X{Et), then for 
any initial condition f G Di^L^), the function u{t,x) = Ex[f{Zt)] solves the higher 
order initial value problem (13.131) . 

Proof. The corollary follows immediately from Theorem 13.91 together with Theorem 
3.1 in HO]. □ 
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Remark 3.11. An easy extension of the argument for Theorem 13.91 shows that, under 
the same conditions, for any n = 2,3,4, .. . both the initial value problem 



(3.16) ^^^ = Y.f^^Lif{x) + Ly{t,xy, u{0,x) = f{x) 

and the fractional Cauchy problem fl2.4p with (3 = 1/n, have the same unique solution 
given by (12.71) with [3 = 1/n. Hence the process Zt = X{Et) is also the stochastic 
solution to this higher order initial value problem. 

This paper has established a connection between processes subordinated to a Brow- 
nian time subordinator, and to an inverse Levy stable subordinator. A similar but 
weaker correspondence can also be established for stable time subordinators. Let 
T{t)f{x) = Ej.[f{X{t))] be the semigroup of a continuous Markov process X{t) and 
let be its generator. Take St a standard symmetric stable Levy process with index 
< a < 2, so that ^[e*^'^^*)] = e"*''"'". Denote by p{t,x) the density of S{t), and let 
Zt = X{\St\). Nane [17] shows that if a = l/m, where / and m are relatively prime, 
then u{t,x) = Ex[f{Zt)] solves 

Ll'-'fix)-Lyit,xy, uiO,x) = fix) 

x=0 

for any bounded measurable function / in the domain of L^, with D'^f bounded and 
Holder continuous for all multi index 7 such that I7I = 21. 

Suppose that X(t) is a self-similar process with Hurst index H, so that X{ct) ~ 
c^X{t) (equality in distribution). Then it is not hard to check that Zt = X{\St\) is 
also self-similar with Zct ~ c^^°'Zt. If 1 < a < 2 and we take Et to be the inverse or 
hitting time process for a stable subordinator with index (3 = I /a, then Ect ~ c'^Et 
and it follows that the process X{Et) is self-similar with the same index as X(|S'f|). 
Corollary 3.1 in [30] shows that Et has moments of all orders, while the mean of {Stl'' 
diverges for p > a. Hence it seems that these two processes are not equivalent. 

Finally, we note that the equivalence established in this paper does not extend 
to strict subdomains of M*^. Consider the Banach space X = L^{M.'^) and the shift 
semigroup [T{t)f] (x) := f{x + t) with generator L^f = -^f and domain V{Lx) = 
{/ G L^(]R^) : /' G L^(]R+)}. In this case, solutions to the higher order initial value 
problem 

(3.17) |-M(t,x) = ^"ff +L,\{t,xy, u{0,x) = f{xy,x>0 

are not unique. To see this, note that if U\,U2 were any two solutions then u = u\ — U2 
would solve 

d 

(3.18) — M(t, x) = L^'^u{t, x); u{0, x) = 0, x > 

L/t' 
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(-^y^'9^<t,x) = -2j2^^pit,x) 



and uniqueness would require that m = is the only solution to 03.181) . However, 



solves (I3.18P as well, so that solutions to (13.17^ are not unique. On the other hand, 
solutions to the corresponding fractional Cauchy problem fl2.4p with /5 = 1/2 are 
unique [S]. Hence the two forms are not equivalent on this domain. It is an in- 
teresting open problem to find governing partial differential equations for Brownian 
time processes on bounded subdomains. For a typical Markov process with genera- 
tor Lx on a bounded domain, one can solve the boundary value problem ii = LxU; 
u = f on the boundary, as the expectation of X{t) where r is the hitting time at the 
boundary. DeBlassie [23] shows that the analogous result does not hold for iterated 
Brownian motion, and it is likely that a resolution of these problems will require a 
novel approach. 
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